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Abstract 

Droplet breakup is a well studied phenomena in Newtonian fluids. One property of this 
behavior is that, independent of initial conditions, the minimum radius exhibits power 
law scaling with the time left to breakup r. Because they have additional structure 
and shear dependent viscosity, liquid crystals pose an interesting complication to such 
studies. Here, we investigate the breakup of a synthetic nematic liquid crystal (N-(4- 
Methoxybenzylidene)-4-butylaniline), known as MBBA. We determine the phase of the 
solution by using a cross polarizer setup in situ with the liquid bridge breakup apparatus. 
Consistent with previous studies of scaling behavior in viscous-inertial fluid breakup, when 
MBBA is in the isotropic phase, the minimum radius decreases as 7-I 03^0.04 nematic 
phase however, we observe very different thinning behavior. Our measurements of the thin- 
ning profile are consistent with two interpretations. In the first interpretation, the breakup 
is universal and consists of two different regimes. The first regime is characterized by a 
symmetric profile with a single minimum whose radius decreases as 7-i-5i±0 06 r^^^ second 
and final regime is characterized by two minima whose radii decrease as 7-0-52±o.ii These 
results are in excellent agreement with previous measurements of breakup in the nematic 
phase of liquid crystal 8CB and 5CB. Interestingly, we find that the entire thinning behav- 
ior can also be fit with an exponential decay such that Rmin ^ exp (^(1.2x102 Hz)Ty This 
dependence is more reminiscent of breakup in polymers where entropic stretching slows the 
thinning process. An analogous mechanism for slowing in liquid crystals could arise from 
the role played by topological constraints governing defect dynamics. Consistent with this 
interpretation, crossed polarizer images indicate that significant alignment of the liquid 
crystal domains occurs during breakup. 
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I. INTRODUCTION 



One well-studied property of Newtonian fluids is their behavior under droplet 
breakup. As a thread of fluid thins until the point of breaking, the minimum thread 
radius decreases to zero in flnite time so that the induced pressure diverges. The 
balance between the gradient of this divergent pressure and terms in the Navier- 
Stokes equation that match this divergence determine the limiting behavior of the 
fluid. The various stress balances that can be acheived in Newtonian fluids have 
been weU studied both in experiments and theory [3 El [JOl [H ESI ESHSU ESI 

ESIEHISISZIIISIISQ]. 

Universality in droplet breakup is not a property of all fluids for any set of envi- 
ronmental and initial conditions. For example, it has been shown that low viscosity 
drops breaking inside a viscous fluid, such as water inside viscous oil [38] or air in- 
side water [28l ESI HH], deviate from universal behavior. In such cases, the initial 
and boundary conditions set before breakup are encoded in the asymptotic behav- 
ior. Deviations from universality have also been studied in non-Newtonian fluids. 
In particular, polymeric systems have been shown to break slower than a power law 
fluid [H |34j- Exponential models seem to flt the data well in these systems where 
elastic stresses eflFectively increase the extensional viscosity [12] . 

More recently scientists have begun to explore how liquid crystal solutions break. 
These systems present an additional challenge since the resulting dynamics arise 
from coarse grained averaging of small scale flows that are intimately related to the 
liquid crystalline structure [H EJ [SI O ESI US]. The details of how these small scale 
structures lead to the bulk flows are poorly understood in part because it has been 
diflicult to image the evolution of the liquid crystalline domains during breakup. 

Recent studies have shown that universal behavior can describe breakup dynamics 
in liquid crystals, but the exact form depends on the phase [HISlESj. In particular, for 
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liquid crystals where the domain size is on the same order as the droplet dimensions, 
it was found that thinning of the minimum droplet radius depends not only on the 
phase but specifically on the orientation of the nematic directors inside the fiuid [U EJ 
[6l [43]. Furthermore, it was shown that this thinning can be described in a universal 
framework in terms of a power law dependance of the minimum radius Rmin as a 
function of time left until breakup r. 

Studies have also been conduced for liquid crystals in which the characteristic 
domain size is much smaller than the drop radius [9l [26]. For example, studies con- 
ducted with the thermotropic liquid crystal 8CB (4-cyano 4-octylbiphenyl) and the 
non-ionic surfactant CuEq (hexaethylene glycol monododecyl ether) in the smectic 
phase, showed universal behavior that could be described by a similarity solution for 
a strain thinning power-law fiuid [32[ |33[ |39l [44] . There, it was shown that the entire 
breakup was symmetric with a single minimum whose radius Rmin oc r°-^. Using 
a cross polarizer setup for the C12EQ studies, it was observed that during breakup 
the domains within the smectic liquid crystal became more aligned and coarsened - 
a well established mechanism for strain thinning. Moreover, it was shown that the 
universal behavior persisted until the neck size was comparable to the liquid crystal 
domain length scale, which could be set by the sample preparation [26j. 

Studies of breakup in 8CB and 5CB (4-cyano 4pentylbiphenl) in the nematic phase 
showed two different breakup regimes [9j. Initially the breakup was symmetric and 
chracterized by Rmin ^ ^^'^ [9^^. As breakup proceeded, the minimum bifurcated, 
and the profiles became asymmetric with Rmin ^ ^^'^ - These exponents also showed 
a slight dependence on temperature. When these liquids were heated so that they 
underwent a transition into the isotropic phase, it was found that breakup proceeded 
similarly to Newtonian fiuids in the inertial- viscous regime [50] where both the initial 
symmetric and final asymmetric profiles are described by Rmin ^ ^^'^ [9j. 

The drawback in the 5CB and 8CB samples is that they are opaque and so the 
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cross polarizer technique of [26] used to to study C^Eq has never been used to study 
nematic hquid crystal breakup. In addition, the thin necks that characterize the 
final nematic regime make collapsing the profiles onto a universal curve prohibitive. 
In the present study, the nematic liquid crystal MBBA (N-(4-Methoxybenzylidene)- 
4-butylaniline) is used to combat the first of these challenges, namely the inability 
to use cross polarizers to probe domain orientation. In addition to studying domain 
orientation, we determine the phase of the solution with the cross polarizer setup in 
situ with the liquid bridge breakup apparatus. We determine the minimum radius 
dependence on phase and analyze the emergent asymptotic behaviors. 

II. SAMPLE 

We use MBBA (98% Sigma-Aldrich) a translucent yellow thermotropic liquid 
crystal with rod-like molecules that have an aspect ratio of about 2.5 [14j (Fig. [T]). 
At room temperature it has a density of 1.03g/cc and a surface tension of about 
30 dyn/cm [37j. Additionally, it has been observed that the surface tension parallel 
to the direction of alignment is greater by about 8 dyn/cm than the surface tension 
perpendicular to the alignment [37] . There are many characterizations of the viscosity 
of a nematic fiuid and a complete description for MBBA has been determined in the 
literature [llj. For example, the Miesowicz coefficient 773, which corresponds to the 
usual shear viscosity of an isotropic ffuid fT4] is about 0.04 Pa • s [16] for MBBA at 
room temperature or roughly 40 times that of water. 




FIG. 1. Chemical structure of MBBA. 



5 




FIG. 2. Phase transition of MBBA under a cross polarizer. Below ^ 20°C, MBBA is a 
solid. In this series of images, the sample begins melted in the Nematic phase (left image) 
and transitions into an isotropic phase (right image). Note that the exact melting point is 
difficult to identify because it is very sensitive to the level of water in the air. 

Under vacuum, at about 20° C, MBBA melts from a crystalline solid into a 
nematic liquid crystal. Then, at about 50° C, the sample melts into an isotropic 
fluid [40j. Unfortunately MBBA is very sensitive to the presence of water and will 
hydrolyze into p-anisaldehyde (pAde) and p — n butylaniline (pBa) [3j, which alters 
its transition temperatures. In fact this property prevented it from being widely used 
in LCD screens [15]. By maintaining a constant temperature and humidity in the lab 
we found that the transition temperatures could be stabilized so that the complete 
phase behavior of MBBA was accessible. 

Since MBBA is translucent it is possible to probe the phase behavior using cross 
polarizers. For example. Fig. [2] shows a sample of MBBA under cross polarizers 
during its phase transition from a nematic to an isotropic liquid at 35° C (Supple- 
mentary Movie 1). Here the nematic portions of the sample appear bright while the 
isotropic portions are opaque. In the final frame on the right, the backlight can not 
be seen as the liquid is now in an isotropic phase. Over all, these in-situ observations 
provide a new handle for studying breakup in nematic liquid crystals since they en- 
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able accurate determination of the phase boundaries associated with the measured 
breakup dynamics under laboratory conditions. 

III. SETUP AND PROCEDURE 

During the experiments a fluid sample is placed between two conducting metallic 
plates that can be separated vertically to induce breakup (Figure |3]). Regulated 
resistors attached to the plates with thermally conductive paste were used to control 
the temperature. In order to minimize heat loss, Non-conducting plates were used to 
separate the conducting plates from the rest of the apparatus. Finally, the resistors 
selected enabled us to generate temperatures between ^ 20° C and 60° C. 

A fast camera was used to image the breakup at a high frame rate of 22099 fps 
with a fleld of view of 128 x 512 pixels, a medium frame rate of 8213 fps with 512x512 
pixels, and a low frame rate of 4796 fps with 800 x 600 pixels. A bellows was used 
to magnify the image at all frame rates so that 1 pixel was about 3 /xm. The sample 
was illuminated with a high luminosity LED light and cross polarizers were used to 
visualize domain orientation as illustrated in flgure [3| 

A Matlab code was used to extract the minimum radius as a function of time 
left to breakup. The various stages of the analysis are shown in Fig. [4] To extract 
the bridge proflle from the individual image frames, we used edge detection software 
(Fig. |4]B) followed by a clustered component analysis that removed noise (Fig. |4]C). 
The extracted proflles were flt to a high order polynomial from which the minimum 
radius is extracted at sub-pixel resolution (Fig. [4]D). To determine r we need to 
identify the time of breakup, which is bounded between two successive frames. To 
mediate this uncertainty, we image many breakup events. By overlaying the points 
from each movie on top of each other and adjusting the time of the flrst frame in 
each movie we have reduced this ambiguity. 
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FIG. 3. Schematic diagram of the apparatus. The sample fluid is suspended between two 
heated plates. The upper plate can be adjusted vertically to induce a narrow fluid bridge 
suitable for breakup. A diffuser is wrapped around the bridge in order to increase the area 
of illumination. Domain orientation in the sample is probed with polarizers on the camera 
lens and touching the diffuser. The bellows facilitate magnification of the point of breakup 
in the sample. 

In addition to overlaying points from movies taken at the same frame rate and res- 
olution, we image the breakup using low, medium and high frame rates to probe the 
small and large r behavior. The data from these regimes are then stitched together 
using the procedure outlined in the previous paragraph to form a single breakup 
curve. To identify scaling behaviors, power law fits to data subsets were made. In 
addition, we performed linear fits to the Rmin versus r data in semi-log to identify 
potential exponential behavior. We found that the largest source of uncertainty in 

8 



FIG. 4. The major steps of the analysis program. Frame A shows the unaltered image 
captured by the camera. Frame B illustrates the data points extracted using the Matlab 
edge identification algorithm. Using a clustering algorithm we are able to remove the 
speckle noise resulting from the edge detection program, shown in frame C. Finally, the 
profile is fit using a high order polynomial and the minimum radius is measured. Shown in 
frame D is the extracted profile overlaid onto the original image along with an illustration 
of the radius we measure. 

the fitting arose from the choice of fit range. To estimate this uncertainty we deter- 
mined how the power law exponents and exponential rates extracted from the fits 
varied as the range of data was decreased by 30%. 

IV. RESULTS 

The dependance of the minimum profile radius for MBBA as a function of time 
left until breakup was examined at 10 different temperatures that span the isotropic 
to nematic phase transition. At each temperature we characterized the dependence 
of the minimum radius on the time left to break up r for 6 movies: three at the 
high frame rate, two at the medium frame rate and one at the low frame rate. For 
example, the minimum radius versus r for nematic MBBA at T = 24.2° C is shown 
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in Figures [5]A.&B. To test for universal behavior, we plot Rmin versus r on a log-log 
plot (Figure [5]A.). We find that the data are consistent with two diflFerent power law 
regimes. Figures [5p IV, V and VI, show three images from the first regime that 
is characterized by symmetric profiles with single minima whose radii decreases as 
^1.6 Figures [5]C II and III show profiles from the second and final regime that is 
characterized by two minima whose radii decrease as r^-^^. The fact that the second 
regime has a lower power law indicates that the dynamics are slowing in the second 
regime. To test whether such dynamics are exponential in form as is the case in 
polymeric systems, we plot Rmin versus r in semi-log form (Figure [5]B). We find that 
both the symmetric and asymmetric regimes can be fit with a single exponential 
form Rmin ^ exp ^(1.3 x 10^ Hz)r^ . Similar behaviors were observed for all 5 bridge 
profile measurements of MBBA in the nematic phase. 

At higher temperatures where MBBA was determined to be isotropic the breakup 
behavior was markedly different than breakup in the nematic regime. For exam- 
ple, the minimum radius versus r for isotropic MBBA at T = 52.4° C is shown 
in Figures |6]A&B. From Figure |6]A, we conclude that in contrast to the nematic 
breakup, the power law behavior consists of a single asymptotic exponent. As seen 
in Figures [6]C II-V for both the initial symmetric ([6]CII, III) and final asymmetric 
([ojCIV, V) regimes the minimum radius decreases as r^-^^ . This exponent and the 
final breakup profiles are consistent with the behavior observed in drop breakup in 
the viscous-inertial asymptotic regime of Newtonian fiuids [50j. We observe that 
unlike for the nematic regime shown in Figure [5]B, the data are inconsistent with an 
exponential dependance of Rmin on r in the isotropic regime (Figure |6]B). 

We summarize the temperature dependence of the power-law exponents and ex- 
ponential rates in Figures [7] and [8| respectively. We find that the power law behavior 
for all ten temperatures in both the nematic and isotropic regimes shows no signif- 
icant temperature dependence (Fig. [7|. Averaging over the different measurements 
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FIG. 5. A&B. Minimum radius i2min as a function of time left until breakup r in the 
nematic phase in a log- log format (A) and semi- log (B). Each plot shows data from six 
movies that have been stitched together. The ambiguity in the r of the first frame was 
reduced by aligning the six curves overlaid on top of each other. The behavior for r > 0.002 
is a transient regime that depends on the initial geometry more than any final scaling. In 
A, we found two power law regimes, one with Rmin ^ ^^'^ and then Rmin ^ ^^'^ in the final 
regime. In B, we observe a single exponential slope of 1.3 x 10^. CI-CVI. Representative 
profiles at T = 24.2° C. While we observe a single minimum in profiles CIV-CVI there are 
two minima in the later profiles CII and CIII. The breakup terminates in a single satellite 
droplet. 
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FIG. 6. A&B. Minimum radius i?min as a function of time left until breakup r in the 
isotropic phase in a log- log format (A) and semi- log (B). Each plot shows data from six 
movies that have been stitched together. The ambiguity in the r of the first frame was 
reduced in aligning the six curves overlaid on top of each other. The behavior for r > 0.002 
is a transient regime that depends on the initial geometry more than any final scaling. In 
A, we found a single power law regime with Rmin ^ t. In B, we observe that an exponential 
model does not well describe the data. CI-CVI. Representative profiles at T = 24.2° C. CII 
and GUI show the initial symmetric behavior while GIV and GIV show the final asymmetric 
breakup. 
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FIG. 7. Shown here are the power laws extracted from the breakup profiles of six movies 
for each temperature. In the nematic regime (gray shaded region), there are two different 
power laws. The larger power laws correspond to the initial symmetric behavior where the 
bridge has a single minimum (e.g. Fig. [5]C IV, V and VI) while the lower branch documents 
the behavior in the final asymmetric regime where the bridge has two minima (e.g. Fig.[5]C 
II and III). Finally, at temperatures above 35° C where birefringence data shows MBBA is 
isotropic, we find a single power law exponent describes both the symmetric and asymmetric 
regimes. 



we find that in the nematic phase, the first breakup regime is characterized by a 
symmetric profile with a single minimum whose radius decreases as t-I-^i^o o^. The 
second and final regime is characterized by two minima whose radii decrease as 
^o.52±o.ii the isotropic phase in both the initial symmetric and final asymmetric 
regimes the minimum radius decreases as 7-i o3±o.o4 pij^^Uy, we note that in both 
the isotropic and nematic regimes, the breakup is followed by the formation of a 
large satellite droplet, seen in Figures [5pl and Figure [6pl. Similarly, the fitted expo- 
nential rates show no significant temperature dependence and are characterized by 



Rrmn ^ exp (^(1.2 ± 0.1 X 10^ Hz)r) (Fig. 



8) 
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FIG. 8. Shown here are the slopes (in Hz) of the exponential models fit to Rmin versus r in 
the semi-log representation, i.e. the coefficient of the linear term in log(i?^^^) = ci log t+cq. 
We find that this fit well-describes the data in the nematic regime (c.f. Fig. [5^) but does 
not model the data well in the isotropic regime (c.f. Fig. [6^) and so the plot is redacted 
to show only the nematic behavior. For each temperature, the slope was extracted from 
the breakup of six movies. 



Since MBBA is translucent, we are able to probe the breakup behavior under cross 
polarizers. Figure [9] shows sample profiles in the nematic regime (room temperature) 
at different times prior to breakup. In Figure |9]A., we find that the structural domains 
in the bridge are randomly oriented. In contrast, in Figure [9]B, we observe large 
scale alignment of the domains to a direction that is parallel with the surface of 
the fluid. This alignment initiates early in the breakup process. Therefore, it will 
be important to take into account the rheology arising from domain realignment in 
order to determine the forcing terms that govern the breakup dynamics. 

Unfortunately, we could not examine the profiles under cross polarizers for small 
r. Given the current camera sensitivity, the main limitation was in the intensity of 
light that illuminated the sample. Since MBBA is translucent and not transparent. 
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FIG. 9. Breakup under cross polarizers. To improve the image, we artificially lightened 
the background (uniform light gray). We note however that the features within the drops 
were not altered. Profile A was imaged at 5.9 seconds prior to breakup. At 0.2 seconds 
before breakup (profile B), one can find that the domains start aligning well before the 
final regime is reached. 



a large backlight intensity is required, so to produce Figure |9] we had to compensate 
by increasing the exposure time at the cost of a lower frame rate. One could try 
to add more LEDs at a higher intensity, but this is nontrivial because order of 
magnitude estimates show that we need to increase the intensity by about a factor of 
100 in order to raise the frame rate enough to see the asymptotic breakup regimes. 
A possible method for achieving this intensity is with a laser. For example, with 
phase-modulated flow birefringence it should be possible to determine the orientation 
angle at a point [23l EU [H]. However, such a method will only report on the local 
orientation. Therefore, obtaining more comprehensive information about domain 
behaviors and alignments would require scanning at multiple positions in the bridge 
and accounting for the deflections of the incident beams as the bridge breaks. 
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V. CONCLUSIONS 



The tabulated results in Figures [7| and |8] confirm that the asymptotic breakup 
behavior of MBBA depends on the liquid crystal phase. In the universal interpre- 
tation, the power law changes as the fluid transitions between the isotropic and 
nematic phases. As was noted in [9j, the flnal power law regime in the nematic phase 
is consistent with the 2/3 power law expected in the inertial Newtonian regime. 
However, qualitative differences in the proflle suggest very different radial and axial 
scalings [42] and therefore rule out this scaling relation. In the nematic phase, the 
asymptotic behavior is described by two power laws which correspond to qualitative 
changes observed in the droplet proflle. In contrast, in the isotropic regime, the 
data is well modeled by a single exponent for both the symmetric and asymmetric 
flnal breakup proflles. Furthermore, this behavior is consistent with observed drop 
breakup in the viscous-inertial asymptotic regime of Newtonian fluids. Bulk rheol- 
ogy measurements indicate Newtonian behavior with a constant viscosity in both 
the isotropic and anisotropic regimes for shear rates of up to 2 kHz [27j. Never- 
theless, strain thinning behavior has been shown to be important for determining 
breakup similarity solutions observed in smectic liquid crystals [26j. One possible 
resolution to these contrasting observations is that the breakup flows may be explor- 
ing substantially higher shear rates where non-Newtonian rheologies might become 
important. 

In the exponential interpretation, the temperature dependance is clear from the 
fact that the asymptotic behavior has a different relationship between Rmin and r 
across the phase transition. The universal scaling in the isotropic regime is consistent 
with expectations. Lacking the additional structure from long range molecular orien- 
tations, MBBA in the isotropic regime resembles a Newtonian fluid in its asymptotic 
behavior. In contrast, the asymptotic behavior of MBBA in the nematic phase is well 
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described by a single exponential curve. Furthermore, this relationship seems to hold 
across various temperatures, as summarized in Fig. [8| The exponential model has a 
clear motivation that can be understood in an analogy to polymeric systems. Studies 
have shown that the power law dependance of Rmin on r in a Newtonian fluid can 
change abruptly to an exponential relationship when polymers are added ^ El |34j . 
This transition is due to the elasticity of the polymer molecules. The increasing 
energy required to stretch the polymers is realized as an increase in the effective ex- 
tensional viscosity [12j. In a similar way, there is an increasing energy cost in liquid 
crystal systems. With nematic ordering, there is a preferred long range order for the 
molecules within a liquid crystal domain. Near breakup, topological constraints may 
require the molecules to rearrange. This rearrangement can lead to elastic stresses 
and slower thinning. 

While the current dataset cannot differentiate between the power law scaling and 
the exponential model, the later has a clear motivation. In the universal interpre- 
tation, the measurements for the the nematic and isotropic regimes are consistent 
with the corresponding scalings found in the same phases of 5CB and 8CB [9, 26j, 
indicating that similar dynamics govern all three systems. Therefore, if further stud- 
ies reveal that exponential slowing is a more appropriate model for the asymptotic 
behavior of these systems, we expect this conclusion to hold more generally. Higher 
precision measurements of the asymptotic breakup behavior coupled with further 
analysis of domain orientation will elucidate the underlying power law or exponen- 
tial dependance. In addition, such studies will be able to pinpoint the mechanism 
that is driving the observed flows. 

In addition to bulk rheology measurements of MBBA, studies have been conducted 
to investigate the average orientation of the nematic directors with respect to the 
flow direction. We have already noted that for nematic fluids with domain sizes on 
the order of the droplet dimensions this orientation of nematic directors influences 
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the breakup dynamics in addition to the thermodynamic phase [H El [6l [43]. In our 
analysis of breakup in MBBA, where the domain size is much smaUer than the initial 
droplet size, we have probed the dependance of scaling only on the phase. There 
are several methods for extending the study to additionally probe the dependance on 
director orientation. For example, there is a known relationship between temperature 
and preferred director angle with respect to the flow direction [131 113 SSI SB] • If 
the measurement uncertainties can be reduced it may be possible to elucidate how 
director alignment alters asymptotic breakup. Additionally, external electric flelds 
have been related to director orientation [27j, and could therefore also be used as a 
control knob for tuning alignment and possibly breakup behavior. Having established 
a clear link between the nematic liquid crystal phase and the measured behavior for 
the evolution of the bridge proflle, it will be important to identify and develop handles 
on the speciflc effects of orientation to understand and ultimately control this type 
of breakup behavior. 

One handle that has been utilized for the first time in this experiment is the in situ 
imaging of a nematic liquid crystal breakup under cross polarizers. The analogous 
measurements for a smectic liquid crystal [26] provided insight into the liquid crystal 
orientation and our measurement has initiated this study for the nematic liquid 
crystal domain orientation in MBBA. Future work in this direction will need to 
find ways for increasing the intensity for illuminating the sample, but ideas already 
exist to begin this research. With the added complexity of nematic liquid crystals, 
traditional techniques must be extended in order to fully understand the breakup 
dynamics. However, with tools like in situ imaging under cross polarizers, we can 
begin to probe the deeper relationship between small scale nematic liquid crystal 
structure and bulk fiows during droplet breakup. 
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